Why the magnetic interactions in Naa;Co02 are 3D 
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The puzzle of 3D magnetic interactions in the structurally 2D layered oxide NaccCo02 is addressed 
using first principles calculations and analysis of the exchange mechanisms. The calculations agree 
with recent neutron results, favoring AFM stacking of FM planes. Superexchange via direct 0-0 
hopping and through intermediate Na sp^ hybrids couples each Co to its nearest and six next- 
nearest interplanar neighbors equally. The individual exchange constants are rather 2D, like the 
lattice itself, but due to multiple c-axis exchange paths, the magnetism becomes effectively 3D. 



The layered transition metal oxide (TMO), Naa;Co02 
is attracting considerable interest because of its un- 
usual magnetic and transport properties, and the 
recently discovered superconductivity in its deriva- 
tive, Nai/3Co02-4/3H20. Experimental signatures of 
triplet pairing^'^'^ have lead to suggestions that spin 
fiuctuations5iSiii2i& mediate the superconductivity, in 
analogy with other layered superconducting TMO's like 
high-Tc cuprates and Sr2Ru04, both of which are close 
to magnetic instabilities, as is NaccCo02Ai& However, 
in contrast to those materials where not only the elec- 
tronic structure, but also the magnetic interactions are 
strongly 2D, recent neutron measurements of the magnon 
dispersion^^'^^ in Naa;~o.8Co02 indicate that the anti- 
ferromagnetic (AFM) exchange between C0O2 planes is 
nearly as strong as the ferromagnetic (FM) in-plane ex- 
change. Since dimensionality plays a key role in un- 
conventional superconductivity, it is important to un- 
derstand the microscopic physics behind the apparent 
3D character of magnetic interactions observed in Refs. 
^9, 20 and to clarify the relationship of the magnetic prop- 
erties of the parent compound to those of the supercon- 
ducting hydrate. We address these issues here. 

Na2;Co02 has an unusual magnetic phase diagram as 
a function of x, which itself poses interesting, unre- 
solved questions. For most x, Naa;Co02 is a metallic 
paramagnet, but in a very small range, around x=0.5, 
a charge-ordered, insulating, and possibly antiferromag- 
netic (AFM) region emerges^^. Surprisingly, the metallic 
states on either side of this region are quite different. For 
X < 0.5 the susceptibility is Pauli paramagnet like with 
weak T dependence, while for x > 0.5 it is Curie- Weiss 
like, suggesting local moments. Finally, for x greater 
than ~ 0.75 a spin density wave condenses, ^^iiiSiiiJ^^ with 
clear antiferromagnetism. Tat ~ 22i^ at x=0.82. While 
this and the negative Weiss constant ^'^ii^Ai^ suggest an- 
tiferromagnetic interactions, experiments also show a 
characteristically ferromagnetic hysteresisii as well as 
predicted^ in-plane ferromagnetic fiuctuationsil. 

These varied data are reconciled by neutron scattering 
experiments, ^^'^^ which attribute the Tm =22K transition 
to an A-type AFM ordering (FM planes stacked antifer- 
romagnetically along the c-axis). By fitting magnon dis- 
persion curves to a linear spin wave model, both groups 
conclude that in-plane and perpendicular magnetic ex- 



change constants are of comparable magnitude, indi- 
cating magnetic isotropy, despite considerable structural 
two-dimensionality. We will show here that this quasi- 
isotropy does not imply comparable magnitudes of the 
nearest neighbor exchange constants in- and out-of-plane, 
but rather unexpectedly large coupling to next nearest 
neighbors across the planes, assisted mainly by Na sp^ 
hybrid orbit als. 

Our first principles local density approximation (LDA) 
calculations were done using the augmented plane 
wave plus local orbital (APW+lo) and linearized aug- 
mented planewave (LAPW) methods implemented in 
two codes^i*2S. The experimental lattice parameters, 
a=2.828A, c=10.94A, and LDA relaxed oxygen height 
2^0=0.0859 were used for all calculations.^^ The partial 
occupation of Na was treated both by the virtual crystal 
approximation (VGA) and in supercells. In the VGA, 
each 2d Wyckoff position of the PQ^/mmc cell was oc- 
cupied by a fictitious atom, atomic number, Z = 10 + 
to model a partial occupancy of x. We also did some 
calculations with Na at the 2h positions. Supercell cal- 
culations were done at selected x with real Na at some 
2d or 2h positions. 

Prior calculations show that there is a FM solution in 
the LDA^^ for Naa;Go02 for all x in the experimentally 
relevant range. Our calculations for a tripled x = 2/3 
supercells show that the A-type AFM order recently re- 
vealed by neutrons is actually the preferred LDA ground 
state. AFM stacking of the ferromagnetic layers is fa- 
vored by 2.3 meV/Go for Na in the 2h site (Na on top 
of Go) in a x supercell, and 1.7 meV/Go in the 
2d site. This is consistent with a recently measured^^i 
metamagnetic AFM to FM transition at the relatively 
low field of 8T. We find very similar results in the VGA: 
1.4 meV/Go with Na in site 2d and 2.2 meV/Go with Na 
in site 2h. The total spin moment inside the Go APW 
spheres for the A-type AFM ordering is only ~ 0.02 jiB 
less than that in the FM case which is half metallic. 
In other words, the Go is maximally spin polarized for 
the doping level x consistent with the ~ 0.2 jiB hmit 
from neutron^^ and muon spin relaxation^^ experiments 
at X ~ =0.8. However, note that at x = 2/3, samples are 
paramagnetic; the LDA ordered ground state is presum- 
ably suppressed by quantum critical fluctuations that are 
potentially important for superconductivity.- 
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We now turn to the relative magnitude of in-plane and 
perpendicular exchange constants. Since the band struc- 
ture is metallic, antiferromagnetic superexchange com- 
petes with ferromagnetic double exchange (the former de- 
pends on hopping linearly, and the latter quadratically). 
The large in-plane dispersion leads to a net FM in-plane 
interaction, while the net inter-plane coupling, due to the 
smaller c-axis dispersion, is AFM. As mentioned, nearly 
isotropic 3D magnetic interactions are unexpected in lay- 
ered compounds. Still this is not inconsistent considering 
the large bonding-antibonding splitting of the aig band 
in LDA calculations^^ at large x. At x=2/3, this split- 
ting at r is 0.21 eV, i.e. 15% of the full t2g bandwidth. 
If the band structure is mapped onto an effective Co- 
only model and only nearest neighbor hopping across the 
planes is allowed, as assumed in Refs. 19 and 20 (we will 
argue that this is not a good approximation), and the 
effective hopping amplitude is t±^ then the aig splitting 
at r would be 12t± and there would be no first order 
splitting for the two bands (consistent with the LDA 
band structure). This gives t± ^ 15 meV. Although the 
t2g — t2g iu-plauc hoppiug is nearly an order of magni- 
tude larger, this still gives less anisotropy than would be 
anticipated for a layered material. 

To proceed, we need to understand the physics of the 
interlayer coupling. The hopping must proceed, via O p 
states, which, in turn, requires either direct O Pz — Pz 
hopping or hopping assisted by diffuse unoccupied Na s 
and p orbitals. In the latter case the details of Na place- 
ment may be important. Moreover, it is significant that 
the energy separation between the unoccupied Na 3s and 
3p states in Naa;Co02 is found to be rather small, smaller 
in fact than the corresponding band width. This allows 
the s and Px^y orbitals of a Na atom, sitting inside an 
Oq prism, to combine and form bands that may be de- 
scribed as coming from sp"^ hybrid orbitals, specifically, 
hi = {s- V2py)/VS, /i2,3 = {s± VV^Px -^Py/V2)/V3. 
These are asymmetric, orthonormal, and directed to the 
midpoints of the O prism edges. So the Na-assisted part 
of the 0-0 hopping goes from one O to another O above it 
via one, two, or three Na sp'^ hybrids depending on how 
many of the three nearby Na sites are occupied. Note 
that if only Na s states were involved, the hopping am- 
plitude to a second nearest neighbor O in the next plane 
via a specific Na atom would be the same as that of hop- 
ping to the O right on top. Thus Na p participation 
changes the interlayer coupling in an essential way. 

Besides the Na assisted hopping, direct 0-0 hopping 
may be important. Considering the rhombohedral O site 
symmetry, the relevant hopping is from Co aig to O Pz- 
Further, in NaajCo02, the O atoms in adjacent C0O2 
sheets are directly on top of each other, favoring inter- 
layer hopping via the opposing O Pz orbitals. To es- 
timate this contribution we suppressed the Na-assisted 
hopping in three ways. First, we looked at C0O2 + x 
e~ i.e., a compound without Na, but with the correct 
number of valence electrons compensated by a uniform 
background charge. At x = 2/3, the resulting aig band- 




FIG. 1: (color online) Charge density of the aig bonding 
orbital in X2/3C0O2 with X= e" (left), X= VGA Na (middle), 
and X= real Na (right). The Na ions facilitate 0-0 bonding 
and thus superexchange. The VGA and supercell calculations 
differ in the asymmetry of the sp^ hybrid. 



splitting at F is 0.12 eV or ^^70% of the total effective 
coupling in Na2/3Co02. We verified that this is inde- 
pendent of the amount of extra valence charge in the 
system, by doing calculations with 0.82, 0.75, and 0.4 
extra e~ . Next, for x = 2/3, we put Ne in the Na sites, 
again adding electrons with a compensating background. 
Finally, we did calculations for Na2/3Co02, where Na 
denotes Na with s and p orbitals artificially shifted to 
higher energy. The results change very little, indicating 
that the splitting from direct 0-0 hopping is ~ 0.12 eV. 
Thus, the effective interplanar coupling can be written as 
t± = t^^{x) + tj, where t^^{x) is a function of doping 
and is constant, and t^^{x)/t^ varies from 0.48 to 
0.78 in the range 0.6 < x < 0.9. 

The Na assisted hopping can be assessed using VCA 
or supercell calculations. These differ in two ways. First, 
in the VCA, all Na sites are occupied regardless of x, 
while in a supercell and in reality a portion, 1 — x are 
empty. The number of effective hopping paths via Na, 
and, correspondingly, the effective hopping between the 
planes, should therefore be reduced in this proportion 
relative to the VCA. Second, the VCA nuclear charge 
is Z = 10 + X, so the unoccupied Na s and p energies 
are moved up. This will reduce the Na-0 hybridization, 
which is inversely proportional to the energy separation 
between Na 3s{p) and O 2p bands. This implies a arti- 
ficial monotonic increase of Na-assisted hopping with x 
in the VCA. These two effects are in opposite directions 
and should at least partially cancel. To assess this, we 
performed VCA calculations for < x < 1, and supercell 
(real Na) calculations for x=0, 1/2, 2/3 and 1. The cal- 
culated VCA aig splitting at F depends on x non-linearly 
in both cases; the supercell calculations give similar re- 
sults with slightly larger splittings. (Fig. (2)). We verified 
that the deviation from this linearity for the supercell cal- 
culations has the same origin as in the VCA: the Na 3s 
and 3p levels shift with x though to a lesser degree than 
in the VCA (the shift is due to the changing Coulomb 
potential as a function of doping). This is seen in the 
charge density of the bonding aig state at F, (Fig. QJ. 
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FIG. 2: A comparison of VGA and supercell aig splitting at 
the r point, with the same structural parameters (see text), 
showing a monotonic increase with x. Galculations without 
Na orbitals, but with the proper valence charge show that the 
0-0 hopping contribution is independent of x. 

The relatively weak O Pz — Pz overlap is augmented in 
the VGA by a composite state in the middle of the O- 
O bond, consisting of the three Na sp^ hybrids; removal 
of one of these with the corresponding Na atom in the 
Na2/3Co02 supercell makes the O orbitals tilt towards 
the remaining Na ions. 

A somewhat counterintuitive consequence of the above 
analysis is the fact that the nearest neighbor approxima- 
tion is invalid for inter-plane hopping. Neither a Go-0- 
Na-O-Go path nor a Go-O-O-Go path needs to end on the 
Go site directly above the one where it started. There are 
9 different Na-containing paths connecting nearest neigh- 
bor Go ions in different planes^^, and 3 connecting one Go 
with each of the 6 second neighbors^^. For paths with- 
out Na, there are 3 that connect nearest neighbors and 1 
that connects second neighbors. Assuming that each Na- 
containing path contributes r to the effective interplanar 
hopping amplitude and each 0-only path contributes r', 
we find that k)/t2(k) = Srx/r' and 

t2(k) =r (3 + 2cosak+2cosbk+2cosck), (1) 

where a, b, and c are respectively the three projections 
onto the ab plane of the vectors connecting a given Go ion 
with the three O above. At F, = 27rx, and t^=0.12 
eV, thus we find, at x=0.82, r = 0.5 meV and = 1.1 
meV. Finally, we use the number of r and r' paths that 
connect a Go ion to each type of neighbor to calculate 
tc and t^, the hopping integrals for first and second Go 
interplanar coupling to get: 

Je _ _ f QT-X + 3T' Y _^ 

We now revisit the interpretation of the spin wave dis- 
persion observed in Refs. and ^20, including second 
neighbor assisted hopping. A straightforward generaliza- 
tion of Eq. 2 in Ref. [jjl yields 



E = 2SJ{M0) - J\\ (0) + J\\ (q) + ( — )}2 - {JM)Y 




FIG. 3: Magnon dispersion of NaccCo02 along high symmetry 
lines shown for x= 0.7 and x=0.82. The results of two experi- 
ments in which data were fit to a nearest-neighbor interaction 
model are compared to our dispersion which includes second- 
nearest out-of-plane neighbors. Our parameters were fit to 
the data of Ref. for x=0.7, adjusted for Na content (see 
text), and used to get the dispersion for x=0.82 without re- 
fitting. 



where, following their notation, hkl is the wave vector in 
units of the reciprocal lattice vector. J^(q) = ^c(q) + 
J"^'(q) and J"^'(q) = 2J^ cos(7r/)[2 cos(27r/i) + 2 cos(27r/c) + 
2 cos(27r(/i+/c)]. Using the k = (0, 0, 1/2) magnon energy, 
we can extract Jc = 1.98 meV and =0.22 meV. This 
fits experiment well with Jy = —4.5 meV. Bayracki et 
al^ using nearest neighbors only, found Jc = 3.3 meV and 
J|| = —4.5 meV, leading to the conclusion of magnetic 
isotropy. Fig. El shows that our model produces results 
that are completely indistinguishable from the nearest 
neighbor model used in Ref. [l^ however, the parameters 
are more physically reasonable and compatible with the 
highly anisotropic electronic structure of the compound. 

Finally, we address the differences in the magnon dis- 
persions obtained in Ref. 19 and Ref. 20, both shown 
in Fig. |31 Because of notational differences between Ref. 
[T^ and Ref. l2dL we extract exchange constants directly 
from the Ref. [23 data, using the same formalism as Ref. 
19 and get Jc = 6.1 eV and Jy = -4.5 eV. Using VGA 
calculations at x=0.7 and x=0.82, we scale the parame- 
ters of our nearest and next-nearest neighbor model and 
compare them with the results of Helme et alm^ We es- 
timate that at X = 0.75 Jc=1.69 meV, =0.19 meV, 
and J|| = —4.5 meV. So, the in-plane hopping remains 
unchanged, in agreement with experiment, and the inter- 
planar hopping decreases by ~ 8%, since less Na is avail- 
able to mediate it. 

The measured dispersion^^ in the z direction is much 
larger at x=0.75 than at x=0.S2. Other effects, not ac- 
counted for here, must therefore be operative. One pos- 
sibility is a different pattern of Na ordering. Indeed, in 
our X = 2/3 supercell calculations placing all the Na at 
the 2b sites instead of the 2d sites yields a ~ 35% increase 
in the interlayer coupling. While this number no doubt 
depends on the exact Na arrangement, it does indicate 
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sensitivity to the ordering of a size sufficient to explain 
the experiment. In comparing with experiment, it should 
be also noted that Co disproportionation, seen by NMR 
studies^^ around x=0.7, would change the spin arrange- 
ment of the lattice and therefore the magnon dispersion. 
Co ions with formal valency (4-x)+ split proportionally 
into non-magnetic Co^+ and S = 1/2 Co^+. The partic- 
ular arrangement is unknown, but the exchange interac- 
tion may be changed. Note that the phase boundary of 
charge order is not yet clearly defined and calculations^^ 
suggest that, near the crossover, there is a region where 
two distinct magnetic ions exist. Another consideration 
is that the spin- wave model adopted here and in Refs. ASk 
and[2nl assumes a rigid spin moment of S=l/2, which may 
not be a good approximation for the weak magnetism in 
this compound. The actual moment, both measured and 
calculated, is smaller and grows with x. Last but not 
least, one should keep in mind that FM and AFM in- 
teractions here are competing with each other thereby 
amplifying the relative effect of doping changes on mag- 
netism. 

In summary, we find that inter-plane coupling to sec- 
ond neighbors plays an important role in the out-of-plane 
magnetism of Na^cCo02. Heisenberg type models includ- 
ing only one nearest neighbor exchange across the plane 
are insufficient; exchange between the six next nearest 
neighbors in adjacent planes is needed. Fortunately, one 
can estimate the ratio of the two exchange constants. 



making it possible to extract them from experiment with- 
out increasing the number of fitting parameters. The re- 
sulting exchange constants provide a magnon spectrum 
that matches experiment extremely well and yields a 
physically realistic picture of magnetic interactions in 
this layered material. The LDA energy difference be- 
tween the FM and the A-type AFM ordering for Na 
content x = 2/3is~2 meV/Co, which is reasonable 
considering the exchange coupling deduced from experi- 
mentally measured spin wave dispersion (6.6 meV and 
12.1 meV/Co in Refs. 19 and 20, respectively). We 
note that superexchange drops with the distance more 
strongly than other coupling mechanisms, such as dou- 
ble exchange. Thus, although the inter-planar coupling 
is surprisingly strong in Naa;Co02, yielding three dimen- 
sional magnetic character, it would be expected to be 
very weak in the hydrated superconducting compound. 
It is even possible that a cross-over from AFM to FM 
coupling occurs with hydration due to suppression of the 
superexchange interaction. If the magnetic interactions 
become effectively 2D in the hydrated compound, fiuc- 
tuations would be enhanced, possibly revealing super- 
conductivity associated with a nearby magnetic quantum 
critical point. 
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